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Twenty years ago Poole et al. (Nature 360, 324, 1992) suggested that the anomalous properties of supercooled water 
may be caused by a critical point that terminates a line of liquid-liquid separation of lower-density and higher-density 
water. Here we present an explicit thermodynamic model based on this hypothesis, which describes all available ex- 
perimental data for supercooled water with better quality and with fewer adjustable parameters than any other model 
suggested so far. Liquid water at low temperatures is viewed as an 'athermal solution' of two molecular structures with 
different entropies and densities. Alternatively to popular models for water, in which the liquid-liquid separation is 
driven by energy, the phase separation in the athermal two-state water is driven by entropy upon increasing the pres- 
sure, while the critical temperature is defined by the 'reaction' equilibrium constant. In particular, the model predicts 
the location of density maxima at the locus of a near-constant fraction (about 0.12) of the lower-density structure. 



Cold and supercooled water have been the subject of inten- 
sive experimental, theoretical and computational studies for 
the last several decades. Still, the famous statement of 1972 
by Franks 'of all known liquids, water is probably the most 
studied and least understood' 1 remains topical. This is espe- 
cially true for metastable supercooled water, which is now a 
focal point of debates. On the other hand, there is a growing 
interest in the prediction of properties of supercooled water. In 
particular, in applied atmospheric science it is commonly ac- 
cepted that the uncertainties in numerical weather prediction 
and climate models are mainly caused by poor understanding 
of properties of water in tropospheric and stratospheric clouds, 
where liquid water can exist in a deeply supercooled state 2 ' 3 . 
Reliable prediction of properties of supercooled water is also 
important for cryobiology 4 . 

A provocative, but thermodynamically consistent, view on 
the global phase behaviour of supercooled water was ex- 
pressed in 1992 by Poole et al. 5 . According to this view, 
deeply in the supercooled region, just below the line of ho- 
mogeneous ice nucleation, there could exist a critical point 
of liquid-liquid coexistence (LLCP) that would terminate 
the line of first-order transitions between two liquid aque- 
ous phases, low-density liquid (LDL) and high-density liquid 
(HDL) (Fig. 1). The anomalies in the heat capacity, the com- 
pressibility, and the thermal-expansion coefficient experimen- 
tally observed upon supercooling 6 " 13 thus might be associated 
with this critical point, even if it is inaccessible. 

Intriguing liquid-liquid phase separation and the existence 
of the second critical point in water still remain a plausible 
hypothesis which needs further verification. In view of the 
inaccessibility of the LLCP to direct experiments, develop- 
ment of an equation of state, based on a solid physical con- 
cept and able to accurately describe all available experimental 
data, might help in resolving the supercooled-water dilemma. 

In this paper, we offer an approach to thermodynamics 
of phase separation in supercooled water alternative to com- 
mon views. Following Bertrand and Anisimov 20 , we assume 
that liquid water is a non-ideal athermal 'solution' of two 
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FIG. 1. Hypothetical phase diagram of cold water. Supercooled 
water exists between the melting temperature 14 Tm and the homo- 
geneous ice nucleation temperature 15 7h- Below 7h, there may ex- 
ist a liquid-liquid critical point, marked by 'C\ which terminates a 
liquid-liquid coexistence curve. The location of this critical point 
is shown as predicted by our model. The adopted location of the 
liquid-liquid coexistence curve is close to similar suggestions of 
Mishima 13 ' 16 and Kanno and Angell 7 . The extension of the coex- 
istence curve into the one-phase region is the line of maximum fluc- 
tuations of the order parameter, the Widom line 17 . Thin solid lines 
represent phase boundaries between the different ices 18 ' 19 . 



supramolecular arrangements, which undergoes phase separa- 
tion driven by non-ideal entropy upon increase of the pressure. 
In the athermal two-state model, the non-ideality driven by 
entropy determines the critical pressure, whereas the critical 
temperature follows from the condition of 'chemical reaction' 
equilibrium. We have developed an explicit equation of state, 
which is based on the athermal two-state model and which 
describes all available experimental data for supercooled wa- 
ter, both H2O and D2O, with better quality and with fewer ad- 
justable parameters than any other model suggested so far. We 
have also shown that the popular two-state 'regular-solution' 
scenario, in which phase separation is energy-driven, fails to 
describe the experimental data on water. 
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I. HOW A PURE LIQUID CAN UNMIX 

The only possible fluid-fluid phase transition in simple one- 
component substances, such as argon or methane, is separa- 
tion into liquid and gas. This transition is driven by attrac- 
tion forces between molecules. The critical temperature of the 
liquid-gas separation is uniquely proportional to the depth of 
the intermolecular interaction potential, which is a superpo- 
sition of the attraction and repulsion potentials. The simplest 
model of the liquid-gas transition is the lattice gas, which, 
despite its simplicity, reflects the most important features of 
fluid phase behaviour 21 . A central concept in physics of phase 
transitions is the notion of the order parameter 22 . The order 
parameter in fluids is the difference between the densities 
of liquid and vapour. Thus, in simple one-component fluids 
only one liquid phase and only one critical point are possible. 
To observe a liquid-liquid separation one needs a mixture of 
two or more species. Remarkably, the same lattice-gas model 
that describes the liquid-gas transition in pure substances can 
also explain liquid-liquid demixing in weakly-compressible 
binary solutions. In binary liquid solutions the order param- 
eter is the difference in concentrations of two liquid phases, 
while the critical temperature is proportional to the difference 
in the depths of the interaction energies between like and un- 
like molecules. 

However, such a simple fluid-phase behaviour is not a law 
of nature. As explained by Mishima and Stanley 23 , if the in- 
termolecular potential of a pure fluid could exhibit two min- 
ima, the interplay between these minima may define the criti- 
cal temperature T c and pressure P c of liquid-liquid separation. 
Application of this picture to water implies that the second, 
liquid-liquid, critical point is driven by molecular interaction 
energy just like the liquid-gas critical point. 

The hypothesized existence of two liquid states in pure wa- 
ter can be globally viewed in the context of polyamorphism, 
a phenomenon that has been experimentally observed or the- 
oretically suggested in silicon, liquid phosphorus, triphenyl 
phosphate, and in some other molecular-network-forming 
substances 24 - 25 . Commonly, polyamorphism in such systems 
is described as energy-driven. However, there is an ambigu- 
ity in terminology adopted in Refs. 24 and 25, where the 
term 'density, entropy-driven' is used for an energy-driven 
phase separation. There is an example of entropy-driven liq- 
uid polyamorphism, isotropic Blue Phase III in chiral liquid 
crystals 26 . 

The thermodynamic relation between the molar volume 
change AV and the latent heat (enthalpy change) of phase tran- 
sition, Q = TAS (where AS is the molar entropy change) is 
given by the Clapeyron equation dP/dT = Q/TAV = AS/AV. 
Therefore, the relation between the volume/density (p = 1 /V ) 
change and the latent heat/entropy change is controlled by the 
slope of the transition line in the P—T plane. 

Two features make the second critical point in water phe- 
nomenologically different from the well-known gas-liquid 
critical point. The negative slope of the liquid-liquid phase 
transition line in the P—T plane in supercooled water (Fig. 1) 
means, in accordance with the Clapeyron equation, that the 
higher-density liquid water is the phase with larger entropy. A 



large value of this slope at the critical point (about 30 times 
greater than that for the gas-liquid transition at the critical 
point) indicates the significance of the entropy change with 
respect to the density change, and, correspondingly, the im- 
portance of the entropy fluctuations. Secondly, supercooled 
water tends to separate upon pressurizing. The relative sig- 
nificance of the entropy change, combined with a high degree 
of cooperativity of hydrogen bonds 27 , suggest that the liquid- 
liquid phase separation in water near the LLCP may be driven 
by entropy rather than by energy. 

The Clapeyron equation itself does not provide an an- 
swer whether the liquid-liquid transition in pure substances 
is energy-driven or entropy-driven. To answer this question 
one should examine the source of non-ideality in the free 
energy. Bertrand and Anisimov 20 introduced the concept of 
a 'lattice liquid', as opposed to a lattice gas, an imaginary 
one-component liquid which exhibits liquid-liquid separation 
upon pressurizing with a vertical liquid-liquid transition line 
in the P—T plane, thus without the density change. The or- 
der parameter in a lattice liquid is the entropy and the 'order- 
ing field', conjugate to the order parameter, is the tempera- 
ture. The phase transition in a lattice liquid is purely entropy- 
driven. The critical temperature (the same as the temperature 
of the liquid-liquid transition line) is defined by thermody- 
namic equilibrium as zero ordering field, but not by the in- 
teraction potential. In real water, as seen in Fig. 1, the slope 
of the hypothesized LDL-HDL transition line in supercooled 
water changes from very steep at higher temperatures to rel- 
atively flat at lower temperatures. It is thus tempting to as- 
sume that the liquid-liquid separation in water may represent 
a special kind of liquid polyamorphism, intermediate between 
two limiting cases: mostly entropy-driven phase separation 
(lattice liquid-like) near the critical point and mostly energy- 
driven (lattice gas-like) separation into two amorphous states 
observed in water at about 140 K 28 . We show in this paper that 
the actual behaviour of supercooled water appears to be much 
closer to the lattice-liquid behaviour than to the lattice-gas be- 
haviour. 

The entropy-driven separation of a lattice liquid can be fur- 
ther specified in terms of a two-state model 20 . The pure liquid 
is assumed to be a mixture of two interconvertible states or 
structures of the same molecules, whose ratio is controlled by 
thermodynamic equilibrium. The existence of two states does 
not necessarily mean that they can separate. If these states 
form an ideal 'solution', the liquid will remain homogeneous 
at any temperature or pressure. However, if the solution is 
non-ideal, a positive excess Gibbs energy, G E — H E — TS E , 
would cause phase separation. If the excess Gibbs energy is 
associated with a heat of mixing H E , the separation is energy- 
driven. If the excess Gibbs energy is associated with an ex- 
cess entropy S E , the separation is entropy-driven. The entropy- 
driven nature of such a separation means that the two states 
would allow more possible statistical configurations, and thus 
higher entropy, if they are unmixed. 

Two-state models for liquid water have a long his- 
tory, dating back to the 19th century 29 - 30 . More recently, 
two-state models have become popular to explain liquid 
polyamorphism 24 ' 25 ' 31 . Ponyatovsky et al? 2 and Moynihan 33 
assumed that water could be considered as a 'regular binary 
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solution' of two states, which implies that the phase separa- 
tion is driven by energy. None of these two-state models for 
water has been used for a quantitative description of available 
experimental data. 



II. VIRTUAL LIQUID-LIQUID CRITICALITY IN 
SUPERCOOLED WATER 

There is no direct experimental evidence of the LLCP in real 
water, but it is indirectly supported by thermodynamic ar- 
guments based on density measurements 13 and by critical- 
like anomalies of thermodynamic response functions 6-8 . The 
known existence of two states of glassy water 28 , as well 
as experiments on nano-confined water, which does not 
crystallize 34 ' 35 , are also consistent with the possibility of 
a 'virtual' liquid-liquid separation. This possibility is also 
supported by simulations of some water-like models, such 
as ST2 36 . The exact location of the liquid-liquid critical 
point in these models is uncertain 27 . Moreover, for the mW 
model 37 it has been recently shown that spontaneous crys- 
tallization occurs before a possible liquid-liquid separation 
could equilibrate 38 ' 39 . However, the anomalies observed in the 
metastable region of the mW model 37 " 40 might still be asso- 
ciated with the existence of a virtual LLCP in the unstable 
region. 

The first attempt to develop an equation of state for super- 
cooled water, based on the assumption that the LLCP exists, 
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and on the asymptotic theory of critical phenomena 
made by Fuentevilla and Anisimov 43 and further elaborated 
and clarified by Bertrand and Anisimov 20 . In particular, both 
works estimated the LLCP critical pressure below 30 MPa, 
much lower than most of simulated water-like models pre- 
dicted. Holten et al. 44 used the same asymptotic equation of 
state, also in a mean-field approximation 45 , but introduced 
the noncritical backgrounds of thermodynamic properties in 
a thermodynamically consistent way. The resulting correla- 
tion represents all available experimental data for supercooled 
water, H2O and D2O, within experimental accuracy, thus es- 
tablishing a benchmark for further developments in this field. 
However, there is a concern regarding the application of the 
asymptotic theory to a broad range of temperatures and pres- 
sures including the region far away from the assumed critical 
point. Such an extension makes the description of experimen- 
tal data inevitably semi-empirical since all non-asymptotic 
physical features are absorbed by the adjustable backgrounds 
of thermodynamic properties. This fact underlines the need 
to develop a closed-form theoretically-based equation of state 
which would satisfy the asymptotic critical anomalies and, at 
the same time, describe regular behaviour far away from the 
critical region. 



III. WATER AS AN ATHERMAL SOLUTION OF TWO 
STATES 

We assume pure liquid water to be a mixture of two intercon- 
vertible states or structures A (HDL) and B (LDL). The frac- 
tion of molecules in state B is denoted by x, and is controlled 



A^B. 



(1) 



The states A and B could correspond to different arrange- 
ments of the hydrogen-bonded network. 46 The Gibbs energy 
per molecule G is the sum of the contributions from both 
states, 



BA 



G = (1 -x)n A +xii B = fi A +xfi 



where fi A and /i B are the chemical potentials of A and B in the 
mixture. The variable x is conjugate to /Z BA = /X B — ji A . If A 
and B form an athermal non-ideal solution, the Gibbs energy 
of the mixture is 

f /-A z^-BA 

\j \j \j 

= \-x hxlnx+ (1 — jc) ln( 1 — x) + axil —x), 

k B T k B T k B T V ' y ' V ; ' 

(2) 

and the chemical-potential difference is then 



M BA = G BA. 



k R T 



ln- 



1 -x 



■<o(l-2x) 



where G BA = G B — G A is the difference in Gibbs energy per 
molecule between pure configurations A and B, k B is Boltz- 
mann's constant, and co = C0(P) is the interaction parameter, 
which depends on pressure but not on temperature. 

Considering x, the fraction of B, as the 'reaction coordinate' 
or 'extent of reaction' 47 , the condition of chemical reaction 
equilibrium, 



dG 



= m ba = o, 



T.P 



yields the equilibrium constant K = K(T,P) of reaction (1) as 

G BA 



InK 



k B T 



and defines the equilibrium fraction x e through 

Xc 



lnA' + ln- 



1 -X e 



■0)(1 -2x e ) =0. 



(3) 



The expression (2) combined with the equilibrium condition 
x — x e is our equation of state for the two-state supercooled 
water. The non-ideality of the two-state mixture is entirely as- 
sociated with the excess entropy of mixing S E = — Q)x(l —x), 
while the heat (enthalpy) of mixing is zero. The parameter 
co determines the critical pressure through its pressure depen- 
dence, which we approximate as 

co = 2 + coqAP, 

where AP = (P — P c )/Pck B T Q with a subscript 'c' denot- 
ing critical parameters. An alternative, regular-solution model 
would have a factor w/k B T in place of CO in equation (2), 
where the interaction parameter w would determine the crit- 
ical temperature. For a regular solution, the excess entropy is 
zero, while the heat of mixing is H E = wx{ \ —x), thus making 
phase separation purely energy driven. 
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The conditions for the critical point of liquid-liquid equi- 
librium, 

( dx l ) T.P ' ( dx l ) T,P 

yield in the case of the athermal-solution model the critical 
parameters x c = 1/2 and G)(P C ) = 2. At pressures below the 
critical pressure, CO < 2 and the Gibbs energy versus fraction 
x — x e shows a single minimum. Above the critical pressure, 
CO > 2 and there are one or two minima in the Gibbs energy, 
depending on the value of the equilibrium constant K(T,P). 
If there are two minima, the minimum with the lowest Gibbs 
energy represents stable equilibrium, and the other minimum 
corresponds to a metastable state. When lnK(T,P) = 0, both 
minima have the same Gibbs energy, representing two-phase 
equilibrium. The critical temperature T c and the location of 
the liquid-liquid transition (LLT) curve in the P—T plane are 
thus determined by the dependence of the equilibrium con- 
stant K on temperature and pressure. Since the LLT curve is 
defined as an analytical function of temperature and pressure 
lnK(T,P) = 0, this function is to be obtained from the exper- 
imental data. In Supplementary Section 1 we explain how we 
match the LLT curve of the two-state model to the experimen- 
tally expected shape and location. The resulting expression for 
InK is 

InK = X (AT + aAP + bATAP), 

where AT = (T — T c )/T c . The parameter a = — p c k^dT /dP 
is the slope of the LLT curve at the critical point, and b de- 
termines the curvature. The parameter X is proportional to 
the heat of reaction (1), while the product V = ka is pro- 
portional to the volume change of the reaction. Since X/v = 
AS/pckftAV = Q/TAV, the heat of reaction is asymptoti- 
cally related to the latent heat of phase separation through 
XAx = Q/k#T c , while Ax = AS/k^X — Q/Xk^T c serves as the 
order parameter along the LLT. 

The Gibbs energy G A of the pure structure A defines the 
background properties and is approximated as 

G A = Y,c mn (Af) m (APT, (4) 

where m (0 to 3) and n (0 to 5) are integers and c m „ are ad- 
justable coefficients. Expressions for thermodynamic proper- 
ties which follow from our model are given in Supplementary 
Section 2. 



IV. EFFECTS OF CRITICAL FLUCTUATIONS 

Thermodynamics predicts the divergence of fluctuations at 
the critical point. Entropy fluctuations are proportional to 
the isobaric heat capacity Cp, volume/density fluctuations are 
proportional to the isothermal compressibility Kj, and cross 
entropy-volume fluctuations are proportional to the thermal 
expansion coefficient (Xp 12 . The two-state model described 
above is essentially mean-field, not affected by fluctuations. 
Being expanded near the critical point in powers of x — x c , 
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FIG. 2. Optimization of the LLCP location. The coloured map 
shows the reduced sum of squared residuals. The solid red line is 
the hypothesized LLT curve. The dashed curve shows the tempera- 
ture of homogeneous ice nucleation 15 . The blue dotted curve is the 
LLT suggestion by Mishima 13 and the green dotted curve is the 'sin- 
gularity' line suggested by Kanno and Angell 7 . 

equation (2) takes the form of a Landau expansion 22 . The pro- 
cedure to include the effects of critical fluctuations is well 
developed and known as crossover theory 42 . This procedure 
is fully described in Supplementary Section 3. The variables 
x — x c and P — P c are renormalized such that the behaviour 
close to the critical point agrees with the asymptotic critical 
behaviour 44 , and crosses over to mean-field behaviour given 
by the equation of state, (2) and (3), far away from the critical 
point. 

V. DESCRIPTION OF EXPERIMENTAL DATA AND 
DISCUSSION 

We have fitted both the mean-field and the crossover formu- 
lation of our two-state model to experimental data, as de- 
scribed in Supplementary Section 4. We adopt the location 
of the liquid-liquid coexistence parallel to the homogeneous 
nucleation curve, which appears to be close to a suggestion of 
Mishima 13 ' 16 based on the shape of metastable melting curves 
of different ices, and a 'singularity line' of Kanno and Angell 7 
based on the extrapolation of the compressibility anomalies. 
The optimum locations of the critical point form a narrow 
band in the P—T diagram. The best fit for the critical point 
is obtained at about 227 K and 13 MPa, with X = 2.3 and 
C0o = 0.35. The LLT curve was chosen to intersect the band 
of LLCP locations at the optimal value of the critical pressure 
(Fig. 2) where it has a dimensionless slope of I /a = 15.3. This 
location of the critical point is about 10% of p c ksT c higher 
than that optimized by the mean-field approximation (see Sup- 
plementary Section 4A). This shift is induced by critical fluc- 
tuations, as follows from crossover theory. 

The description of the density and the response functions 



5 



1200 




200 220 240 260 280 300 
Temperature (K) 



FIG. 3. Density of cold and supercooled water as a function of 
temperature and pressure. Black curves are the predictions of the 
crossover two-state model. Tyi (thin black) indicates the melting tem- 
perature and Tn indicates the homogeneous nucleation temperature. 
The thin blue line is the liquid-liquid equilibrium curve, with the 
critical point C. The red line is the line of maximum density, and the 
green line is the line of a constant LDL fraction of about 0.12. Sym- 
bols represent experimental data 12 ' 13,48 . Mishima's data 13 have been 
shifted by at most 0.3% to bring them into agreement with data for 
stable water, as described in Ref. 44. 



are shown in Figs. 3 and 4. A significant improvement com- 
pared to the previous works 20 ' 44 ' 45 is that our equation of state 
does not show any sign of additional thermodynamic insta- 
bility beyond the liquid-liquid separation. We now believe 
that the mechanical and thermal instability below the LLT, 
reported previously 20 - 44 ' 45 , were associated with the asymp- 
totic nature of the equation of state, which caused nega- 
tive backgrounds of the compressibility. The number of ad- 
justable background coefficients in equation (4) is now 14, to 
be compared with 16 44 ' 45 (see Supplementary Table SI). An- 
other improvement is the agreement of the coexistence densi- 
ties, shown in Fig. 5, with the experimental densities of low- 
density and high-density amorphous water. We also fitted our 
model to all available experimental data for supercooled D2O 
with the same quality as for H2O (see Supplementary Sec- 
tion 4B). With the same number of adjustable parameters, a 
regular-solution two-state model fails (about ten times higher 



root-mean-square error) to describe experimental data on su- 
percooled water. 

In Fig. 6, we compare the LDL fraction x predicted by our 
equation of state and that obtained in simulations of the mW 
model by Moore and Molinero 40 . Remarkably, both results 
show a similar temperature dependence. In contrast, this frac- 
tion for an ideal LDL/HDL solution would be almost a linear 
function of temperature. The experimental maximum density 
line is located approximately along the line of constant LDL 
fraction x = 0.12 as shown in Fig. 3. 

We confirm the previous finding 20 ' 43 of a critical pressure 
that is much lower than found in simulations 36 . We believe 
that the low critical pressure reflects the entropy-driven nature 
of liquid-liquid criticality in supercooled water. As shown by 
Stokely et alP ', the LLCP critical pressure is determined by 
the ratio of hydrogen-bond cooperativity and hydrogen-bond 
covalent strength. The higher the ratio, the lower the critical 
pressure. 

The necessity to develop a microscopic model for water, 
which would be consistent with its athermal two-state char- 
acter, and which would clarify the microscopic nature of the 
order parameter, (x — x c )/x c , is evident. 

Finally, the problem of the influence of the hypothesized 
liquid-liquid separation on homogeneous ice nucleation re- 
mains unresolved. The fact that the LLT curve is located 
just below the homogeneous nucleation curve and imitates 
its shape suggests that homogeneous ice nucleation may be 
caused by the entropy/structure fluctuations associated with 
the liquid-liquid transition. A connection between the change 
in the structure of liquid water and the crystallization rate of 
ice was shown for the mW model 39 . Whether the LLCP of 
simulated water-like models is in the metastable region or just 
projected to be in the unstable region, the liquid-liquid criti- 
cality still may be responsible for the observed anomalies in 
the accessible domain. 



REFERENCES 



Franks, F. Introduction-Water, the unique chemical. In Franks, F. (ed.) The 

Physics and Physical Chemistry of Water, vol. 1 of Water: A Comprehensive 

Treatise, chap. 1 (Plenum, 1972). 
2 Rosenfeld, D. & Woodley, W. L. Deep convective clouds with sustained 

supercooled liquid water down to -37.5 °C. Nature 405, 440^142 (2000). 
3 Heymsfield, A. J. & Miloshevich, L. M. Homogeneous ice nucleation and 

supercooled liquid water in orographic wave clouds. J. Atmos. Sci. 2335- 

2353 (1993). 

4 Song, Y., Sharp, R., Lu, F. & Hassan, M. The future potential of cryopreser- 
vation for assisted reproduction. Cryohiology 60, S60-S65 (2010). 

5 Poole, P. H, Sciortino, F, Essmann, U. & Stanley, H. E. Phase behaviour 
of metastable water. Nature (London) 360, 324-328 (1992). 

6 Angell, C. A., Oguni, M. & Sichina, W. J. Heat capacity of water at ex- 
tremes of supercooling and superheating. J. Phys. Chem. 86, 998-1002 
(1982). 

7 Kanno, H. & Angell, C. A. Water: Anomalous compressibilities to 1.9 kbar 
and correlation with supercooling limits. J. Chem. Phys. 70, 4008^1016 
(1979). 

8 Kanno, H. & Angell, C. A. Volumetric and derived thermal characteristics 
of liquid D2O at low temperatures and high pressures. /. Chem. Phys. 73, 
1940-1947 (1980). 

'Tombari, E., Ferrari, C. & Salvetti, G. Heat capacity anomaly in a large 
sample of supercooled water. Chem. Phys. Lett. 300, 749-751 (1999). 



6 




200 250 300 200 250 300 200 250 300 

Temperature (K) Temperature (K) Temperature (K) 

FIG. 4. Response functions as a function of temperature, a, Isothermal compressibility 713,49 . Pressures, from top to bottom: 0.1 MPa, 
10 MPa, 50 MPa, 100 MPa, 150 MPa, 190 MPa, 300 MPa, 400 MPa. b, Heat capacity at constant pressure Cp (open circles: data from Archer 
and Carter 10 , closed circles: data from Angell et al. 6 ) and at constant volume Cy (calculated) at 0.1 MPa. See Supplementary Section 5 for 
more details, c, Thermal expansivity 11 ' 12,50 . Pressures, from top to bottom: 380 MPa, 200 MPa, 100 MPa, 0.1 MPa. In a, b, and c, the curves 
are the prediction of the crossover two-state model, and the symbols represent experimental data. 
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FIG. 5. Density along the liquid-liquid transition curve. The 

dashed line represents the mean-field two-state model; the solid 
line represents the crossover two-state model. The open circles are 
the densities of the low-density amorphous (LDA) and high-density 
amorphous (HDA) phases of water at 200 MPa 28 . One can notice that 
the crossover LLT curve is flatter than the LLT curve in the mean- 
field approximation and that the actual position of the critical point 
is shifted to a lower temperature by critical fluctuations. 
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FIG. 6. Fraction of molecules in the low-density state. The fraction 
X is shown for the two-state model at 0. 1 MPa, in the case of an ather- 
mal solution (solid) and an ideal solution (dotted). The dashed curve 
is the fraction of four-coordinated water molecules, i.e., the fraction 
of water molecules with four neighbours, for mW water simulations 
performed by Moore and Molinero 40 (see Supplementary Section 6). 
The temperature is scaled by the temperature of maximum density 
Tmd (250 K for the mW model and 277 K for real water). 
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